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We carry out extensive computer simulations to study the Lyapunov instability of a two- 
dimensional hard disk system in a rectangular box with periodic boundary conditions. The system 
is large enough to allow the formation of Lyapunov modes parallel to the x axis of the box. The 
Oseledec splitting into covariant subspaces of the tangent space is considered by computing the full 
set of covariant perturbation vectors co-moving with the flow in tangent-space. These vectors are 
shown to be transversal, but generally not orthogonal to each other. Only the angle between co- 
variant vectors associated with immediate adjacent Lyapunov exponents in the Lyapunov spectrum 
may become small, but the probability of this angle to vanish approaches zero. The stable and 
unstable manifolds are transverse to each other and the system is hyperbolic. 
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I. INTRODUCTION 



Lyapunov exponents measure the exponential growth, or decay, of infinitesimal phase space perturbations of a 
chaotic dynamical system. For a D-dimcnsional phase space, there are D exponents, which, if ordered according to 

^*\ • size, Xi > Ai+i, are referred to as the Lyapunov spectrum. The classical algorithm for the computation is based 
on the fact that almost all volume elements of dimension d < D in tangent space (with the exception of elements 
of measure zero) asymptotically evolve with an exponential rate which is equal to the sum of the first d Lyapunov 
exponents. Such a d-dimensional subspace may be spanned by d orthonormal vectors, which may be constructed by 
the Gram-Schmidt procedure and, therefore, are referred to as Gram-Schmidt (GS) vectors. The GS-vectors are not 
covariant, which means that at any point in phase space they are not mapped by the linearized dynamics into the 
GS vectors at the forward images of that point [l|. As a consequence, they are not invariant with respect to the 
time-reversed dynamics. Due to the periodic re-orthonormalization of the GS vectors only the radial dynamics is 
exploited for the computation of the exponents, whereas the angular information is discarded. 

Although the angular dynamics is not a universal property and may depend, for example, on the choice of the 
coordinate system [2|, it would be advantageous for many applications, to span the subspaces mentioned above by 
covariant vectors and to study also the angular dynamics of and between these vectors. It has the additional advantage 
to preserve the time-reversal symmetry for these tangent vectors, a property not displayed by the GS vectors. Recently, 
an efficient numerical procedure was developed by Ginelli et al. [l| for the computation of covariant Lyapunov vectors. 
Here we apply their algorithm to a two-dimensional system of rigid disks. 

The choice of hard elastic particles is motivated by the fact that their dynamics is comparatively simple, and 
• • ■ their ergodic, structural and dynamical properties are well known and are thought to be typical of more realistic 

. £h ' physical systems [3j . Secondly, hard-particle systems in two and three dimensions serve as reference systems for 
the most successful perturbation theories of dense gases and liquids [J, lf| . Finally, the combination of a Lyapunov 
analysis with novel statistical methods for rare events Q seems particularly promising for the study of such rare 

5t i transformations in systems, for which hard core interactions are at the root. 

The paper is organized as follows. After an introduction of the basic concepts for the dynamics of phase space 
perturbations in Section UH we summarize in Section IIIII the features and our numerical implementation of the 
algorithm of Ginelli et al. [lj for the computation of covariant vectors and covariant subspaces. In Section IIV1 the 
Henon map serves as a simple two-dimensional illustration. The hard-disk model is introduced in Section [V] In this 
work we restrict ourselves to 198 disks, a number which is dictated by computational economy, but still large enough to 
allow the study of Lyapunov modes. In Section [V Al we study the relative orientations of Gram-Schmidt and covariant 
vectors, which give rise to the same Lyapunov exponents. Next, in Section TVB I we compare the localization properties 
in physical space for these two sets of perturbation vectors. The configuration and momentum space projections of 
the perturbation vectors - Gram-Schmidt or covariant - are the topic of Section IV CI The central manifold (or null 
subspace) and its dependence on the intrinsic continuous symmetries - translation invariance with respect to time 
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and space - is discussed in Section IV Dl Although the null subspace is completely orthogonal to the unstable and 
stable subspaces, it is essential for a proper understanding of the Lyapunov modes 0, @1- Section FVEl is devoted to a 
discussion of these modes and how they are represented by the covariant vectors. In Subsection IV Fl we compute the 
angles between the covariant modes and test for tangency between covariant Oseledec subspaces. In Section IVII we 
conclude with a summary. 

II. PHASE SPACE AND TANGENT SPACE DYNAMICS 

The dynamics of a system of hard disks is that of free flight, interrupted by elastic binary collisions. If r denotes 
the state of the system at time 0, the state at time t is given by T t — <fi f (To), where 0* : X — > X defines the flow in 
the phase space X. Similarly, if STq is a vector in tangent space TX at To, at time t it becomes <5r t = D0*|r o • 5Tq, 
where Dip* defines the tangent flow. It is represented by a D x D matrix, where D is the dimension of phase space. 
A subspace E^ 1 ' of the phase space is said to be covariant if 

Dcf, t \r E^(T Q )=E^(cj ) t (T )). (1) 

This definition also applies to covariant vectors, if E^ 1 ' is one-dimensional. Loosely speaking, covariant subspaces 
(vectors) are co-moving (co-rotating in particular) with the tangent flow. An analogous relation holds for the time- 
reversed flow. 

Next we consider the decomposition of the tangent space into subspaces according to the multiplicative ergodic 
theorem of Oseledec pj|4ll| . Here, we closely follow Ref . Q . 

The first part of the multiplicative ergodic theorem asserts that the real and symmetric matrices 

A ± = lim ([^|r] T ^|r) 1/2 '" (2) 

exist for (almost all) phase points T. Here, T denotes transposition. The eigenvalues of A + are ordered according 
to exp(A' 1 - 1 ) > ••• > exp(A' '), where the A'- 7 ' are the Lyapunov exponents, which appear with multiplicity mP'. 
For symplectic systems as in our case, X^' = — A" +1-J ', which is referred to as conjugate pairing. Similarly, the 
eigenvalues of A_ are exp(— A^') > • • • > exp(— X 1 - 1 '). The eigenspaces of A± associated with exp(±A^) are denoted 
by U± . They are pairwise orthogonal but not covariant. If the A"' are degenerate with multiplicity m^' — dim U± , 
all multiplicities sum to D, the dimension of the phase space. Since the matrices A± are symmetrical, each of the two 
sets of eigenspaces, {U± }, completely span the tangent space, 

TX(r) = U { ± ] (r) © • • • © U ( l ] (T). (3) 

The eigenspaces U± are not covariant, but the subspaces 

U ( l ] © • • • © Uf and U W © • • • © U U \ j € {1, . . . ,£} , (4) 

are. They are, respectively, the most stable subspace of dimension Y7i=j mW °f ^+> an d the most unstable subspace 
of dimension Yh=i m(l) °f ^— (corresponding to the most stable subspace of that dimension in the past). 

The second part of Oseledec' theorem asserts that for (almost) every phase-space point T there exists another 
decomposition of the tangent space into covariant subspaces E^ J '(T) referred to as Oseledec splitting, 

TX(r) = £ (1) (r)©---©£; (£) (r). (5) 

For ST <E E^'(T) the respective Lyapunov exponent follows from 

lim ilog||D#| r -ST\\ = ±X^ Vj€{l,. ..,£}. (6) 

t->±oo |£| 

The subspaces E^> are covariant (see Eq. (p}) but, in general, not orthogonal. According to Ruelle [l(j], they are 
related to the eigenspaces U± of A±: 



E^ = ( t/ W © • • • © U U) 



) n (uf © • • • © uf) . (7) 



This equation is at the heart of the construction of covariant vectors according to Ginelli et al. as described in the 
next section. Furthermore, one can show that 

F U) = E W . . . e E U) = xjM © • • • © U U) (8) 

are covariant subspaces. 



III. NUMERICAL CONSIDERATIONS 

Numerical methods probe the tangent space by a set of D tangent vectors, such that the Lyapunov exponents are 
repeated with multiplicities, Xi > •• • > Xjj. Here, the lower index is referred to as the Lyapunov index. The relation 
between the X^' and A^ is given by 

A W = Ay-(j-l) +1 = • • ■ = Xf(j) , 

where f^> = rrv- 1 ' + • • • + m^' is the sum of all subspace dimensions up to j. 

For notational convenience in the following, the vectors g 3 n , j — 1,...,D spanning the tangent space at time t n , 
are arranged as column vectors of a D x D matrix G n = {gW ■ ■ ■ \gn)- The same convention is used below for other 
spanning vector sets such as G n = (g^\ . . . \g®)a.nd V„ = (v* | . . . \vP\ 

In the classical algorithm of Benettin et al. |12j and Shimada et al. jl 31 ] for the computation of Lyapunov exponents, 
an orthonormal set of tangent vectors G„_i at time t„_i is evolved to a time t„ = t„_i + t, (t > 0), 

where 3„_i is the Jacobian of the evolution map taking the phase space point r„_i at time t„_i to r„ at time t n . 
The column vectors of G n at time t„ generally are not orthonormal any more and need to be re-orthonormalized with 
a Gram-Schmidt procedure. This gives the matrix G n with column vectors {g J } n , which form the next orthonormal 
Gram-Schmidt (GS) basis at time t„. These vectors are pairwise orthogonal but not covariant. Each GS renormal- 
ization step is equivalent to a so-called QR decomposition of the matrix G n , G n — G n R n , where the matrix R n is 
upper triangular |14J . The diagonal elements of R n are required for the accumulative computation of the Lyapunov 
exponents. This procedure is iterated until convergence for the Lyapunov exponents is obtained. 

For the computation of a covariant set of vectors {v J } spanning the tangent space for the phase point T = T(0) 
at, say, time to, Ginelli et al. [l| start with a well- relaxed set of GS vectors at to and follow the dynamics forward for 
a sufficiently long time up to t u = to + lot, storing G n and G n (or, equivalently, R n ) for t n = to + tit, n = 0, • • • , ui 
along the way. At t^ a set of unit tangent vectors {f-'} is constructed according to 

vleSi = s P& n{gl,...,gl} \f j e {1, . . . ,D} , (9) 

which serve as starting vectors for a backward iteration from t^ to time to. The vector v J n will stay in S 3 n at 
any intermediate time t n , because S 3 n is the most stable subspace of dimension j for the time-reversed iteration. 
Arranging these vectors again as column vectors of a matrix V n and expressing them in the GS basis at time t n , one 
has V n = G n C n , where the matrix C n is again upper triangular with elements [CJ; = g l n ■ v J n . If, at any step 

n, C n -\ is constructed from C n according to C n -i — [R n ] C n , Ginelli et al. have shown that V n — J n _iV"„_i 
and, hence, the respective column vectors of this matrix follow the natural tangent space dynamics without re- 
orthogonalization. They are covariant but not orthogonal in general. At this stage of the algorithm, renormalization 
of v J n _ 1 is still required to escape the exponential divergence of the vector norms without affecting their orientation. 

After reaching to at the end of the iteration, the vectors v J point into their proper orientations in tangent space 

such that, according to Eq. ([7|), span [Vq,- ■■ ,Vq ) = .E^'fTo) • • • ® E^ J '(To) is the most-unstable subspace of 

dimension f(i' = m^ 1 ' + • • • + rrv- 1 ' of the tangent space at the space point r , going forward in time. If there are 
degeneracies (as in the presence of Lyapunov modes to be discussed below), the Oseledec subspace E^ J ' is spanned 
according to 

E^=vf U - 1) + 1 ®---®vf (3 \ (10) 

where, as in the following, we omit the arguments for the phase-space point. If there are no degeneracies, v? 3 — W- 3 '. 
Similarly, the Gram-Schmidt vectors may be expressed in terms of the eigenspaces of A_, 

[/ W W ~ 1)+1 ©---©9 /0) . 

For nondegenerate subspaces one finds U_ = gf ' |7I. Ha. fl6l| . 

The drawback of this algorithm for many-particle systems is the large storage requirement for the matrices G n and 
G n (or, equivalently, R n ) for the intermediate times t„ = to + tit, n = 0, • • • ,ui, because r must not be chosen too 
large (containing not more than, say, 20 particle collisions). At the expense of computer time, this can be bypassed 
by storing the matrices only for times separated by, say, lOOr intervals and recomputing the forward dynamics in 
between when required during the time-reversed iteration. In this case, also the phase-space trajectory needs to be 
stored. 



IV. A SIMPLE EXAMPLE: THE HENON MAP 
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FIG. 1: The Henon attractor (black line) and a finite-length approximation of its stable manifold (dotted line) are shown. The 
red vectors are the covariant vectors at the phase point as explained in the main text. The blue vectors are Gram-Schmidt 
vectors. 

To illustrate the foregoing algorithm, we apply it to a simple two-dimensional example, the Henon map [171 ] . 



X n +1 

Vn+l 



a- x n + by n , 



with a = 1.4 and b = 0.3. In Fig. Q]thc Henon attractor is shown (black line), which is known to coincide with its 
unstable manifold. An approximation of the stable manifold is shown by the dotted lines. At the point the initial 
GS basis is indicated by the two orthogonal vectors in blue, where one, as required, points into the direction of the 
unstable manifold. If these vectors are evolved forward in time with the GS method for a few hundred steps, the 
two orthogonal GS vectors at the point ui are obtained. Taking these vectors as the initial vectors v^ and v^, the 
consecutive backward iteration yields the covariant vectors at point indicated in red. As expected, one is parallel 
to the unstable manifold, the other parallel to the stable manifold at that point. 



SYSTEMS OF HARD DISKS 



Now we turn to the study of a two-dimensional system of hard disks in a box with periodic boundaries, where the 
particles suffer elastic hard collisions (without roughness), and move along straight lines in between collisions. The 
case of rough hard disks is the topic of a forthcoming publication [18[ . 

The Lyapunov instability of hard disk systems has been studied in detail in the past |19l - [23| . Here we are mainly 
concerned with the differences encountered with the GS and covariant vectors, which, as we have seen, give rise to 
identical Lyapunov spectra. To facilitate comparison with our previous work, we consider reduced units for which 
the particle diameter a, the particle mass m and the kinetic energy per particle, K/N ', are unity. Here, K is 
the total energy, which is purely kinetic, and N denotes the number of particles. Lyapunov exponents are given 
in units of \J ' K/Nma 2 . If not otherwise stated, our standard system consists of N = 198 particles at a density 



p = N/(L x L y ) = 0.7 and a simulation box with an aspect ratio L y /L x — 2/11, which is periodic in x and y. The 
choice of such a small aspect ratio facilitates the observation of the Lyapunov modes to be discussed later. As usual, 
the total momentum is set to zero. 

The state of the system is given by the coordinates and momenta of all the particles, 

T = {<?„,£»„; n=l,-" ,N}. 

Similarly, an arbitrary tangent vector ST - either a Gram-Schmidt vector g or a covariant vector v - consists of the 
respective coordinate and momentum perturbations, 



ST = {Sq n ,S Pn ; n = !,-■■ ,N}. 



(11) 



The time evolution of these vectors and the construction of the map from one Gram-Schmidt step to the next has 
been discussed before 119, Y2A . 
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FIG. 2: Lyapunov spectrum for the 198 disk system described in the main text. The spectrum calculated in the forward 
direction with the GS method is shown by the blue line, the one calculated in the backward direction with the covariant vectors 
by the red line. Reduced indices i/AN are used on the abscissa. Although the spectrum is defined only for integer i, solid lines 
are drawn for clarity. 

Fig. [5] shows the Lyapunov spectrum for this system computed both in forward direction with the GS vectors 
(blue line) and in backward direction with the covariant vectors (red line). The time of the simulation in the forward 
direction is for to + t u = 2.5 x 10 5 r, where t = 0.6 is the largest interval between two successive Gram-Schmidt rc- 
orthonormalizations, which does not affect the spectrum. The backward simulation is for a time t u — to = 2.5 x 10 4 r . 
The time to (usually of the order of 1 x 10 4 r) is required for the preparation of the relaxed initial state at to- It can 
be observed in the figure that the unstable directions in the future correspond well to the stable directions in the 
past and vice versa. Of course, if the sequence of covariant vectors is followed in the forward direction of time, the 
spectrum is identical to the classical GS results (blue line in Fig. [5]). 



Covariant versus Gram-Schmidt vectors 



Whereas the time evolution of the GS vectors is determined by the exponential growth of infinitesimal volume 
elements belonging to subspaces g 1 © • • • © g l for i G {1, . . . , D} according to exp it J2)=i <\? ) i the growth of an 
infinitesimal perturbation representing a covariant vector v l is directly proportional to exp (t Aj), for all i. Thus, it is 
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FIG. 3: Plot of the scalar product norm, (|g l • v % \), for GS and covariant vectors giving rise to the same Lyapunov exponent 
\i, as a function of i. The line is a time average as discribed in the main text. Left panel: full range of Lyapunov exponents; 
Right panel: enlargement of the central part. 



interesting to compare the relative orientations of respective vectors giving rise to the same exponent. In the left panel 
of Fig. |3] the difference in orientation of the two types of vectors is demonstrated by a plot of \g l ■ v l \ as a function 
of i. The black line is an average over 100 frames separated by time intervals of 250t. Since for tangent vectors only 
their direction and not the sense of direction is important, an absolute value is taken (here and for analogous cases 
below), otherwise the scalar product might average to zero over long times, with equal numbers of vectors pointing 
into opposite directions. For the unstable directions in the left half of the left panel, one observes a rapid decrease 
of the scalar product with i and, hence a rapid increase of the angle between respective covariant and GS vectors. 
This decrease is repeated for the stable directions in the right half of the figure. These two parts are separated by 
the mode region, an enlargement of which is shown in the right panel of Fig. [3] and which will be dealt with in more 
detail below. 

In Fig. 2] we show similar projections (time averages of absolute values of scalar products as before) for selected 
covariant vectors with the whole Gram-Schmidt vector set. One observes that the covariant vectors u J ' belong to the 
GS subspace g 1 ® • • • ® g J , for all j € {1, . . . , 47V} and, thus, give rise to the upper-triangular property of the matrix 
R in the QR-decomposition mentioned above. The curves in the figure strongly depend on the choice of j: 

• If it belongs to the unstable subspace E u but does not represent a Lyapunov mode (top-left panel for j = 370), 
there is no obvious orientational correlation with any of the GS vectors with index i < j. For i = j = 1 corresponding 
to the maximum exponent, the covariant and GS vectors are identical. If, however, the covariant vector represents 
a Lyapunov mode as in the bottom-left panel for j = 393, then its angle with the respective GS vector may become 
smaller, giving rise to a scalar product closer to unity. 

• If the covariant vector belongs to the stable subspace E s but does not represent a mode as for j = 700 in the 
top-right panel of Fig. @J it has non- vanishing components in the GS basis for all i < j with the exception of the zero 
subspace 2N — 2 < i < 2N + 3, which is strictly orthogonal. With the exception of the step at the conjugate index 
i = 4N + 1 — j '■ = 93, the origin of which is not fully understood, there is no indication of orientational correlations 
between the covariant vector with any of the GS vectors for i < j. If, however, the covariant vector represents a 
mode as for j = 400 in the lower-right panel of the figure, there is strong orientational correlation not only with the 
respective GS vector with i = 400, but also with its conjugate pair at 4iV + 1 — i (= 393 in our example). 

It is interesting to note that the leading GS and covariant vectors in the null subspace are always identical (up to 
an irrelevant sign): v 



2N-2 



J2.N-2 



Localization 



The maximum (minimum) Lyapunov exponent is the rate constant for the fastest growth (decay) of a phase-space 
perturbation and is dominated by the fastest dynamical events, a locally-enhanced collision frequency. It is not too 
surprising that the associated tangent vector components are significantly different from zero for only a few strongly- 
interacting particles at any instant of time. Thus, the respective perturbations are strongly localized in physical space. 
This property persists in the thermodynamic limit such that the fraction of tangent- vector components contributing to 
the generation of Ai follows a power law oc N~ v , r\ > 0, and converges to zero for N — > oo [2ll [25H27| . The localization 
becomes gradually worse for larger indices i > 1, until it ceases to exist and (almost) all particles collectively contribute 
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FIG. 4: Time averaged absolute value of the scalar product of selected covariant vectors v 3 (as indicated by the labels) with 
the whole set of Gram-Schmidt vectors g l as a function of i. 



to the coherent Lyapunov modes to be discussed below. Similar observations for spatially extended systems have been 
made by various authors [22|, [23, l28l - l30| , which were consequently explained in terms of simple models [3l|, |32f ■ We 
also mention Ref. |33j |. where the tangent-space dynamics of the first Lyapunov vector g 1 for various one-dimensional 
Hamiltonian lattices is compared to that for the Kardar-Parisi-Zhang model of spatio-temporal chaos. The unexpected 
differences found for the scaling properties are traced back to the existence of long-range correlations, both in space 
and time, in the Hamiltonian chains, the origin of which, however, could not be fully disclosed. The same correlations 
are conjectured to be responsible for a slow 1/vN convergence of Ai towards its thermodynamic limit [33J, which is 
also observed for hard-disk systems [19| . 

Up to now, all considerations concerning localization were based on the Gram-Schmidt vectors. Here, we demon- 
strate the same property for the covariant vectors. According to Eq. dill) we define the contribution of an individual 
disk n to a particular perturbation vector as the square of the projection of ST onto the subspace pertaining to this 
disk, 

Mn = (Sq r f + (Sp n ) 2 - 

Since ST is either a GS vector or a covariant vector both of which are normalized, one has Y] ■, \x n = 1, and /j, n may 
be interpreted as a kind of action probability of particle n contributing to the perturbation in question. It should 
be noted that for the definition of /i„ the Euclidean norm is used and that all localization measures depend on this 
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FIG. 5: Localization spectra W for the complete set of Gram-Schmidt vectors (blue) and covariant vectors (red). The details 
of the hard-disk system are given in Section [V] Reduced indices i/AN are used on the abscissa. In the inset a magnification of 
the central mode region is shown. 



choice. Qualitatively, this is still sufficient to demonstrate localization. From all the localization measures introduced 
[251 [30| , the most common is due to Taniguchi and Morriss [22|, [23[ , 



W=^exp[S], S = 



N 



^ Mn l n Mn 



Here, S is the Shannon entropy for the "probability" distribution /j, n , and (• • ■ ) denotes a time average. W is bounded 
according to 1/JV < W < 1, where the lower and upper bounds apply to complete localization and delocalization, 
respectively. In Fig. [5j we compare W obtained for the full set of Gram-Schmidt vectors (blue curve) to that of all 
the covariant vectors (red curve). The spectra are obtained by identifying ST with all vectors of the respective sets, 
i = 1, • • • ,47V. Not too surprisingly, the localization is stronger for the covariant vectors, whose direction in tangent 
space is solely determined by the tangent flow and is not affected by renormalization constraints. Another interesting 
feature is the symmetry W% = W^M+i-i, which is a direct consequence of the symplectic nature of the flow [34J ]. 



C. Tangent space projections 



It is interesting to see how much the coordinate and momentum subspaces contribute to a particular tangent vector 
ST (see Eq. [TTj) . which may be a Gram-Schmidt vector g l or a covariant vector v l , both associated with the same 
Lyapunov exponent Aj . The time-averaged squared projections of ST onto the coordinate and momentum subspaces 
Q and P, respectively, are given by 



N 



'/</ 



£*£ 



Vp 




(12) 



and are plotted in Fig. [S] for the whole set of Gram-Schmidt vectors, and in Fig. [7J for the whole set of covariant 
vectors, i = 1, • • • , AN. One notes that for the Gram-Schmidt case the contributions of r\ q and rj p to a vector g l and 



its conjugate g 



4N+l-i 



are interchanged, whereas for the covariant vectors v l and v 



4N+l-i 



they are the same. This is 



particularly noticeable for the expanded central regions in the respective right panels of Figs. [5] and [7] 
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FIG. 6: Mean squared projections for the full Gram-Schmidt vector set, i — 1, • • • , 42V, onto the coordinate subspace Q, r/q 
(green line), and the momentum subspace P, jjF (black line), for the 198-particle system defined above. Left panel: full 
spectrum; Right panel: magnification of the central mode-carrying region. 
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FIG. 7: Mean squared projections for the full covariant vector set, i — 1, • • • ,42V, onto the coordinate subspace Q, rf^ v (green 
line), and the momentum subspace P, r/p OV (black line), for the f98-disk system defined above. Left panel: full spectrum; Right 
panel: magnification of the central mode-carrying region. 

D. Central manifold and vanishing exponents 

The dynamics of a closed particle system such as ours is strongly affected by the inherent continuous symme- 
tries, which leave the Lagrangian and, hence, the equations of motion invariant. The symmetries relevant for our 
two-dimensional system with periodic boundaries are the homogeneity of time (or invariance with respect to time 
translation), and the homogeneity of space (or invariance with respect to space translations in two independent di- 
rections) . Each of these symmetries is associated with two vector fields with sub-exponential growth (or decay) and, 
therefore, gives rise to two vanishing Lyapunov exponents 35]. At any phase-space point T, the six vectors span a 
six-dimensional subspace Af{T) of the tangent space TX(T), which is referred to as null space or central manifold. 
This subspace is covariant. If the 4iV components of the state vector are arranged as 



yixi ^y ' 



N N . 1 1 N N\ 

: *lx 5 ^y ' Px ' Py ) ' ■ ' i Px ' Py ) > 



(13) 



10 

the six orthogonal spanning vectors, which are the generators of the elementary symmetry transformations, are given 
by@,IH 



ei = 

e 2 = 

e 3 = 

e 4 = 

e 5 = 

ee = 



(pi,pJ,...,*C>P?5 0,0,...,0,0) 

(1,0,...,1,0;0,0,...,0,0) , 
(0,1,...,0,1;0,0,...,0,0) , 



(0,0,..., 0,0; pl,pl,...,Px,Py) , 



(0,0,. ..,0,0; 1,0,. ..,1,0) , 
(0,0,...,0,0;0,1,...,0,1) . 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



ei corresponds to a change of the time origin, e 4 to a change of energy, e 2 and e^ to an (infinitesimal) uniform 
translation of the origin in the x and y directions, respectively, and e§ and e§ to a perturbation of the total momentum 
in the x and y directions, respectively. The six vanishing Lyapunov exponents are located in the center of the Lyapunov 
spectrum with indices 2N — 2 < i < 2N + 3. The first three of these vectors have non- vanishing components only for 
the position perturbations in the 2iV-dimensional configuration subspace Q, the remaining only for the momentum 
perturbations in the 2iV-dimensional momentum subspace P. They are related by e^ = J e^+3 for k G {1, 2, 3}, where 
J is the symplectic (skew-symmetric) matrix. 

Let us consider the projection matrices a and (3 of the GS and covariant vectors, respectively, onto the natural 
basis, 



a-i,k 



= 9 l e k ; ft,* = »* ■ e fc , fee {!,..., 6} i e {2N - 2, • • • , 2N + 3}. 



For i £ {2N — 2, ... , 2N + 3} these components vanish. Without loss of generality, we consider in the following 
example a system with only N — 4 particles in a periodic box, which is relaxed for t r — 10 6 time units, followed by a 
forward and backward iteration lasting for t u — to = 10 5 time units. Very special initial conditions for the backward 
iteration v]^ = g]^ for i = 1, •• ■ , 4N(= 16) are used. The projections at the time to are given in Table Q] for the GS 
vectors, in Table [TTI for the covariant vectors. 



TABLE I: Instantaneous projection matrix a of Gram- Schmidt vectors (for i G {2N — 2, ■ ■ • , 2N + 3}) onto the natural basis 
{efc, 1 < k < 6} of the central manifold. The system contains N — 4 disks. The powers of 10 are given in square brackets. 



i 


Qi,l 


Cti,2 


Cti/i 


Cti,4 


Qi.5 


Cti,6 


2N-2 


-0.766 


0.582 


0.273 


-0.766[-6] 


0.582[-6] 


0.273[-6] 


2N - 1 


0.256 


-0.114 


0.960 


0.256[-6] 


-0.114[-6] 


0.960[-6] 


2N 


0.590 


0.805 


-0.062 


0.590[-6] 


0.805[-6] 


-0.062[-6] 


2N + 1 


-0.611[-6] 


0.782[-6] 


-0.121[-6] 


0.611 


-0.782 


0.121 


2N + 2 


0.575[-6] 


0.544[-6] 


0.611[-6] 


-0.575 


-0.544 


-0.611 


27V + 3 


-0.543[-6] 


-0.304[-6] 


0.783[-6] 


0.543 


0.304 


-0.783 



A comparison of the two tables reveals the following: 

• The six orthogonal GS vectors g l ; i = 2N — 2, ■ ■ • , 2A^ + 3 completely span the null subspace (the squared elements 
for each rows add up to unity in Table U). The same is true for the six non-orthogonal covariant vectors v z ;i — 
2N - 2, • • • , 2N + 3 in Table M 

• The first three covariant and Gram-Schmidt vectors completely agree. This is a consequence of the special initial 
conditions for the former at the time t u as mentioned above. During the backward iteration the three covariant vectors 
stay in their respective subspaces and remain parallel to the GS vectors (which were stored during the forward phase of 
the algorithm). At t they are still identical to their GS counterparts. The first vectors always agree, v^ ~ 2 — g^ ~ 2 , 
if less special initial conditions conforming to Eq. |H] are used. 



11 



TABLE II: Instantaneous projection matrix matrix /? for the the six central covariant vectors onto the natural basis 
{ejc, 1 < k < 6} of the central manifold. The system contains of TV = 4 particles. The powers of 10 are given in square 
brackets. 



i 


Pi,i 


^ 


ft,a 


A,4 


ft,6 


Pifi 


27V -2 


-0.766 


0.582 


0.273 


-0.766 [-6] 


0.582[-6] 


0.273[-6] 


27V -1 


0.256 


-0.114 


0.960 


0.256[-6] 


-0.114[-6] 


0.960[-6] 


27V 


0.590 


0.805 


-0.062 


0.590[-6] 


0.805[-6] 


-0.062[-6] 


27V + 1 


-0.611 


0.782 


-0.121 


0.611 [-5] 


-0.782 [-5] 


0.121 [-5] 


27V + 2 


0.575 


0.544 


0.611 


-0.575 [-5] 


-0.544[-5] 


-0.611 [-5] 


27V + 3 


-0.543 


-0.304 


0.783 


0.543 [-5] 


0.304 [-5] 


-0.783 [-5] 



• Equivalent components have the same mantissa but may differ by a factors of 10 5 or 10 6 , which are related to the 
duration of the relaxation phase t r and of the forward-backward iteration time t u — to . 

The explanation for this behavior [34| is obtained by a rep eated explicit application of the linearized maps for the 
free streaming and consecutive collision of particles [19J, llj] to the six basis vectors e*. One finds that 



D<j) 



,3+3 



(To) 
(To) 



e J (r t ) , 
t e j (r t ) 



J+3 



(Tt) 



(20) 

(21) 



for j <G {1,2,3}. Eq. (|21l) implies that any perturbation vector with non-vanishing components parallel to e4, es, 
or e§ will rotate towards e\, e2, and e^, respectively. It follows i) that the null subspace Af(T) is covariant; ii) that 
the subspaces A/i = span{ei},A/2 = span{e2J and A3 = span{e3J are separately covariant (from Eq. (f20]>): that, 
as was already noted in Ref. [3], A/" can be further decomposed into the three two-dimensional covariant subspaces 
j\f p = spanjei, e^}, J\f x — span{e2, es}, and j\f y = span{e3, e$}. 



E. Lyapunov modes 



We have seen in Section IV Bl that the perturbation vectors are less and less localized, the smaller the Lyapunov 
exponents become, until they are coherently spread out over the physical space and form periodic spatial patterns with 
a well-defined wave vector k. This collective patterns are referred to as Lyapunov modes. The modes were observed 
for hard particle systems in one, two and three dimensions [3 [20, l22l l23l |36[, for hard planar dumbbells [25|, [33, [38( 
and for one and two-dimensional soft particles [23, |39, [4fJ . A formal classification of the modes is given in Ref. [3] ■ 
Physically, they are interpreted as periodic modulations with wave number (k 7^ 0) of the null modes associated with 
the elementary continuous symmetries and conservation laws. Since this modulation involves the breaking of such 
symmetries, the modes have been interpreteted as Goldstone modes |8J]. Theoretical approaches are based on random 
matrix theory [4j], [42[ , periodic orbit expansion (43[ , and kinetic theory 0, [4J, [45| • 

So far the numerical work on Lyapunov modes has been exclusively concerned with the orthonormal Gram-Schmidt 
vectors {g 4 }, i = 1, • • • , 47V. The purpose of this section is to point out some differences one encounters, if the modes 
for the Gram-Schmidt and covariant vectors are compared. 

Fig. [8] shows an enlargement of the mode-carrying region for the Lyapunov spectrum of Fig. [2] In order to 
emphasize the conjugate pairing symmetry Aj = — A 4 jv+i-i for symplectic systems, conjugate exponent pairs are 
plotted with the same index i on the abscissa, where now i € {1, • • • , 27V}. The open circles are computed from the 
Gram-Schmidt vectors in the forward direction of time, the dots from the covariant vectors during the time-reversed 
iteration. Considering the size of the system (TV = 198), the agreement is excellent. 

The steps in the spectrum due to degenerate exponents is a clear indication for the presence of Lyapunov modes. 
According to the classification in our previous work |7|, the steps with a two- fold degeneracy are transverse (T) 
modes - T(1,0), T(2,0) and T(3,0) from right to left in Fig. [21 Similarly, the steps with a four-fold degeneracy of 
the exponents are longitudinal-momentum (LP) modes - LP(1,0), LP(2,0) and LP(3,0) again from the right. The 



arguments (n x 



account for the number of periods of the sinusoidal perturbations in the x and y directions. Since 



our simulation cell is rather narrow, only wave vectors k parallel to the x axis of the (periodic) cell appear, leaving 
for the second argument [7j. As usual, "transverse" and "longitudinal" refer to the spatial polarization with respect 
to k of the wave-like pattern. 

One of our early observations, which greatly facilitates the classification of the modes for the Gram-Schmid vectors 
[3], is that in the limit TV — > 00 the cosine of the angle O between the 27V-dimensional vectors of the position 
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FIG. 8: Enlargement of the mode regime for the Lyapunov spectrum depicted in Fig. [2] The open symbols indicate exponents 
computed from the Gram-Schmidt vectors, the full dots are for exponents obtained from the covariant vectors. 

TABLE III: Basis vectors of (n x ,0) modes for a hard disk system in a rectangular box with periodic boundaries. We use the 
notation c x = cos(k x x), and s x — sin(k x x), where the wave vector is given by k = (k x , k y ) = (2nn x /L x , 0). Here n x £ {1, 2, 3}. 



n 


Basis of T(n) 


Basis of L(n) 


Basis of P(n) 


(») 


(:)■(:) 


(o)'(o) 


( Px )s x J P *)c x 
\Py J \Py J 



perturbations and momentum perturbations converges to +1 for the smallest positive, and to -1 for the smallest 
negative exponents. See the blue line in Fig. [5J Furthermore, the relation 



dp = C±Sq 



(22) 



holds with known constants C±. This means that these vectors are nearly parallel or anti-parallel for large N and 
that the mode classification may be based solely on 5q. Somewhat surprisingly, this property does not strictly hold 
anymore for the covariant vectors. This is shown by the red line in Fig. [51 where cos(0) is seen to differ significantly 
from ±1 for all i outside of the null subspace (for which 394 < i < 399). Unfortunately, this has dire consequences 
for the representation of the covariant vector modes, since they cannot be purely understood as a vector field of the 
position perturbations only as in the GS case. For the purpose of this paper, however, we restrict to the GS-based 
classification of Ref. [7] . 

Transverse modes are two-dimensional subspaces (for the periodic boundary conditions and a rectangular box), for 
which two orthogonal basis vectors are given in Table IIIII As an example, we show in the panels on the left-hand side 
of Fig. [TOl snapshots of the mode T(1,0) for the index i = 393, namely plots of Sq y as a function of q x (top left), and 
of 5p y as a function of q x (bottom left). The respective plots for the x components fluctuate around zero, as expected, 
and are not shown. Analogous plots for the mode T(2,0) with i = 387 are shown in the panels on the right-hand side. 
The blue points are for GS vectors, the red squares for the respective covariant vectors. It is interesting to note that 
the scatter of the points for the position perturbations is smaller for the covariant modes (red squares) than for the 
GS modes (blue dots). A fit shows that the residuals for the covariant modes are smaller by about a factor of two 
in comparison to Gram-Schmidt. Quite the opposite is true for the momentum perturbations in the bottom row of 
panels. Although the proportionality of Eq. (|2"2")l still holds, the scatter of the red squares for the covariant vectors is 
larger than that of the blue dots for the GS vectors. Such a behavior is always observed and is not simple numerical 
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FIG. 9: Time averaged value of cos(0) = (Sq ■ 5p)/(\5q\\5p\) as a function of the Lyapunov index i for a system with N = 198 
hard disks. Here, Sq £ Q and Sp € P are the 2iV-dimensional vectors of all position perturbations respective all momentum 
perturbation for the Gram-Schmidt vectors g i (blue line) and the covariant vectors Vi (red line). The insets are magnifications 
of the mode-carrying region. 

noise. The reason for this behavior is related to the previous discussion in connection with Fig. [9] and needs further 
clarification. 

Longitudinal (L) and associated momentum (P) modes share the same degenerate Lyapunov exponent A^, and 
generally appear superimposed in experimental vectors. With a rectangular box and periodic boundaries, they form 
four-dimensional LP perturbations. The superposition varies periodically with time. This "dynamics" has been 
identified as a rotation of the pure L and P vectors in the standard frame. For details we refer to previous work in 
Rcf. pj. The patterns for the pure L mode are easily recognizable as sine and cosine functions, but those for the P 
modes are not. As is evident from the spanning vectors for L(1,0) and P(1,0) also listed in Table ITlTI the P modes 
are proportional to the instantaneous velocities of all particles, which does not at all constitute a smooth vector field. 
For a pattern to be recognizable, these velocities need to be "divided out". A full mode reconstruction is required 
as is described for the case of Gram-Schmidt vectors in Ref [3|- Here we carry out an analogous reconstruction in 
terms of the covariant vectors and compare them to the GS modes. In Fig. ITTltwo of the reconstructed patterns for 
L and P modes belonging to the four-dimensional LP(1,0) subspace with indices i <E {388,389,390,391} are shown. 
The blue dots are for GS modes, the red squares for covariant modes. To judge the quality of the reconstruction, we 
have included in the top-right panel also the 5q y -veTsns-q x curve, which vanishes nicely as required. 

For comparison, Fig. [T^] gives results for a completely analogous reconstruction, where instead of the position 
perturbations as in Fig. 111! the corresponding momentum perturbations are used. For this example cosine patterns 
were selected, whereas in Fig. [TT] sine patterns were used. As before, blue dots refer to GS vectors, red squares to 
covariant vectors. 



F. Transversality 



From the Lyapunov spectrum of Fig. [5] and the magnification of its central part in Fig. [5J the following inequalities 
are read off, 



Ai > • • • > A2JV-3 > 



A< c 



> A 2 at+4 > • • • > X4N , 



(23) 



where the equal sign applies for the degenerate exponents belonging to modes. [A ] = is sixfold degenerate in our 
case. Conjugate pairing assures that A^ = — A2jv+i-j- The Oseledec splitting provides us with the following structure 
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FIG. 10: Instantaneous transverse Lyapunov modes T(1,0) for index i = 393 (left panels), and T(2,0) for index i = 387 
(right panels) for the 198-disk system. In the panels at the top (bottom) the {/-coordinate perturbations 5q y ({/-momentum 
perturbations Sp y ) of all particles are plotted as a function of their x coordinate, q x , in the simulation cell. The wave vector is 
parallel to the x axis. The blue dots are for Gram-Schmidt vectors, the red squares for covariant vectors. 



of the tangent space, 



where E u = v 1 



TX = E u © M®E S 



(24) 



,,2Af-3 



and.E s 



.2JV+4, 



v are the covariant stable and unstable subspaces, respectively, 



and Af is the null subspace or central manifold. The question arises whether the system is hyperbolic, which implies 
that the angles between the stable manifold E s and the unstable manifold E M are bounded away from zero for all 
phase points (Due to the existence of a central manifold this is referred to as partial hyperbolicity in the mathematical 
literature 46]). Even more, we may ask whether the angles between all Oseledec subspaces and, hence between all 
covariant vectors, are bounded away from zero for all phase space points. To find an answer to that question, we 
compute in the following the scalar products for all covariant vector pairs and present representative results. (This 
procedure reminds us of the so-called coherence angles introduced by d ' Alessandro and Tenenbaum [43, [48[ , measuring 
the angular distance between a physically interesting direction and the direction of maximum perturbation expansion). 
The lines in Figure 1X31 depict the product norms (\v 3 ■ v l \) for selected covariant vectors v^ with all other covariant 
vectors v z , i ^ j. As before, a time average is performed. The panels on the left-hand side provide three examples for 
v^ from the unstable manifold outside of the mode regime (j = 1,200, and 370 from top-left to bottom-left, respec- 
tively), and similarly on the right-hand side from the stable manifold outside of the mode regime (j = 420, 600, and 
792 from bottom-right to top-right, respectively). One immediately observes that the stable and unstable subspaces 
are not orthogonal. As has been mentioned in Section fVDI and is also convincingly demonstrated in the following Fig. 
[T41 the null subspace Af is orthogonal to both E u and E s . For two covariant vectors from the same subspace, E u or 
E s , however, the scalar product does not vanish indicating considerable nonorthogonality. But at the same time it is 
also well bounded away from unity, which means that the two vectors do not become parallel either. However, one 
possible exception may be the covariant vector pairs (v J , v^ +l ) for adjacent Lyapunov exponents in the spectrum. In 
these cases, the scalar product reaches a pronounced maximum in all of the six panels of Fig. [T3]which may still allow 
these vectors to become parallel occasionally. This will be discussed further below. 
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FIG. 11: Reconstructed position perturbations of the pure L(1,0) mode (top panels) and P(1,0) mode (bottom panels). Only 
the patterns proportional to sm(2-Kq x /L x ) are shown. Blue dots: Gram-Schmidt vectors; Red squares: covariant vectors. For 
details we refer to the main text. 



So far we have considered only vectors v J outside of the mode regime. The case of covariant vectors representing 
modes is treated separately in Fig. [TJ] where, as before, time-averaged scalar product norms (\v^ ■ v l \) for i ^ j 
are plotted as a function of i. The standard deviation is too small to be included in the plots. The panels on the 
left-hand side are for j belonging to unstable transversal modes, the panels on the right-hand side for j belonging 
to unstable LP pairs. The curves for the conjugate stable modes just look like the mirror images around the central 
index. Each vector representing a T or LP-mode has significant contributions to the scalar product only for covariant 
vectors belonging to the same degenerate exponent and - to a lesser extent - the corresponding conjugate (negative) 
exponents (where the latter is not true anymore for the LP(3,0) modes in the bottom-right panel of Fig. 1141 where 
no peak around i = 414 is discernible). 

The covariant vectors belonging to transverse (or to LP) modes span covariant Oseledec subspaces E^> with a 
dimension my' equal to 2 (respective 4). To ease the notation, we refer to them as E^ ' in the following, where X is 
cither T(n x ,0) or LP(n x ,0) with n x <G {1,2,3}. The conjugate Oseledec subspaces, E^ '*, have the same dimension 
and are spanned by the respective conjugate covariant vectors. Fig. [14] shows that the covariant vectors spanning 
any of the subspaces E^ ' or E^ '* have a rather small but finite angular distance and, thus, are transversal. The 
Oseledec subspaces representing modes are themselves transversal to all other subspaces of the Oseledec splitting, but 
to a varying degree. The angular distances in tangent space are generally large except between conjugate subspaces 
E^ ' and E^ '* , for which the scalar products of their spanning vectors may become surprisingly large. 

To check more carefully for transversality even in this case, we show in Fig. [15] the probability distribution for 
the minimum angle between the conjugate subspaces E^ ' and E^ '*, where X stands for the T and LP modes 
as indicated by the labels. This angle $ is computed from the smallest principal angle between the two subspaces 
|49l |50| . If the covariant vectors belonging to E^ ' and E^ '* are arranged as the column vectors of matrices V 
and V*, respectively, the QR decompositions V — QR and V* = Q*R* of the latter provide matrices Q and Q*, 
with which the matrix M = Q T Q* is constructed. The singular values of M are equal to the cosines of the principal 
angles, of which $ is the minimum angle. Since $ is never very small, this method works well and does not need 
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FIG. 12: Reconstructed momentum perturbations for the pure L(1,0) mode (top panels) and P(1,0) mode (bottom panels) 
depicted already in Fig. 1111 Only the patterns proportional to cos(2irq x /L x ) are shown. Blue dots: Gram-Schmidt vectors; 
Red squares: covariant vectors. 



more complicated refinements [494511]. It is seen that all distributions are well bounded away from zero indicating 
trans versality for the respective sub-spaces. 

Finally, we concentrate on the minimum angle between the full unstable subspace E u = v 1 © • • 



conjugate stable counterpart E s = v 



400 



• © v 393 and its 
v 7i>li , using the same method as before. These subspaces include the 
mode-carrying vectors studied before. The probability distribution for the minimum angle is denoted by >F and is 
also shown in Fig. [15] (red line). Also this distribution is well bounded away from zero and indicates trans versality 
between E s and E u . We conclude that for finite N the hard-disk systems are (partially) hyperbolic in phase space. 

In Fig. [T3]it was observed that the scalar products between covariant vectors with adjacent indices are rather large 
and possibly may allow tangencies. To study this point more carefully, we follow a suggestion of G. Morriss and 
consider the angle O = cos -1 \v l ■ v^\ between the vectors v l S F^ ' and v-* e F' ', for which i — j is a specified 
positive integer. The probability distributions for angles with i — j = 1,2,4,20,50 are shown in Fig. [TBI an d for 
i — j = 100, 200, 300 in the inset of the same figure. Whereas the probabilities for i — j > 1 are bounded away from 
zero, the distribution for i — j = 1 seems to converge to zero for O — > 0. 

An even more demanding test is given in Fig 1 171 where the minimum of for given i — j > is plotted as a function 
of i — j. The inset provides a magnification of the most interesting region. One observes that the minimum of the 
angle between covariant vectors specifying Oseledec subspaces with i — j = 1 may indeed become very small, but 
this happens with extremely small probability. Our numerical evidence is consistent with the assumption that the 
angle becomes zero with vanishing probability. 
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FIG. 13: The lines are time averaged (100 frames separated by 150 time units) absolute values of the scalar product of a 
specified covariant vector v 3 with all the remaining covariant vectors Vi-tj as a function of i. Left panels from top to bottom: 
j = 1, 200, and 370 belonging to the unstable manifold; Right panels from bottom to top: j = 420, 600, and 792 from the stable 
manifold. 
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VI. CONCLUSION 

A comparison of the covariant vectors with corresponding orthonormal Gram-Schmidt vectors reveal similarities, but 
also significant differences. The vectors associated with the maximum Lyapunov exponent are identical, v 1 = g , and 
also the leading vectors in the central manifold agree, v 2N ~ 2 = g 2N ~ 2 . All the other corresponding vectors generally 
point into different tangent space directions. Whereas the GS vectors are pairwise orthogonal by construction, the 
covariant vectors are not. Most notably, the perturbation contributions from the particles' positions and momenta are 
significantly different and even exhibit a different symmetry between vectors from the stable and unstable manifold as 
in Fig. [5] For the covariant vectors these contributions agree in accordance with the time-reversal symmetry required 
for them, whereas for the Gram-Schmidt vectors these contributions are interchanged. Another significant difference 
is the degree of localization in physical space for the non-degenerate perturbations. As Fig. [5] shows, the covariant 
vectors are much more localized than the GS vectors in accordance with the fact that they are not dynamically 
constrained by re-orthogonalization. 

From a theoretical point of view, an interesting result is that no tangencies occur between the respective unstable 
and stable manifolds E u and E s . In Fig. [13] the probability distribution ^ of the minimum angle between stable and 
unstable subspaces (including the Lyapunov modes) is well bounded away from zero, and even more so for the vectors 
belonging to unstable respective stable modes. Thus, a hard disk system with N — 198 particles as in our case is 
(partial) hyperbolic for all points in phase space. We even find that all Oseledec subspaces are pairwise transversal 
with non-vanishing angles between them. 

We speculate that for TV — > oo the distribution for the minimum angle between E u and E s may possibly reach the 
origin in Fig. 1151 To clarify this point further studies are required 34]. 

The concept of hyperbolicity is closely linked with the notion of dominated Oseledec splitting for all phase space 
points [46j. We may rewrite Eq. © for the Lyapunov exponents, expressed in terms of the covariant vectors, according 
to 

1 ^1 
Xe = J im m E - ln ll^ T |r(t„) v e (T(t n ))\\, (25) 

n=0 

where t n = nr, and r is the short time interval between consecutive re-normalizations of the covariant vectors. Here, 
Xi is expressed as a time average of a quantity 

A, cov (r(<„)) = im^lr^) v e (T(t n ^))\\, (26) 

which is referred to as local or (time-dependent) Lyapunov exponent, and is a function of the instantaneous phase 
point T(t). The Oseledec splitting is said to be dominated, if the local Lyapunov exponents, when averaged over a 
finite time A, do not change their order in the spectrum for any A l arg er than some finite Ao > [3J]. This is a 
very strong condition on the fluctuations of the local exponents [52], l53f . For symplectic systems it is known that 
the domination of the splitting implies that the system is (partially) hyperbolic [46J] . But it is not clear whether the 
converse is true in our case. The discussion of this point is deferred to a forthcoming publication [34|. 

The number and the dimension of the Oseledec subspaces are constant in phase space. There is no entanglement 
of subspaces, which has been identified as one of the main reasons for the occurrence of well established Lyapunov 
modes [52] ■ We refer to Ref. (53[ for a discussion of a simple but physically-relevant model, for which the dimensions 
of the stable and unstable manifolds frequently change along the trajectory. 

An interesting extension of this work is the study of rough hard particles allowing for energy exchange between 
translational and rotational degrees of freedom [34j. Arguably, this is the simplest model of a molecular fluid. No 
Lyapunov modes are found in this case [54] . An analysis in terms of covariant vectors is presently under way and will 
be published separately. 
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